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Differential measurements of the elliptic (112) and hexadecapole (vi) Fourier flow coefficients are 
reported for charged hadrons as a function of transverse momentum [pr) and collision centrality 
or number of participant nucleons (A^part) for Au+Au collisions at ^/sNN = 200 GeV. The V2,4 
measurements at pseudorapidity \ri\ < 0.35, obtained with four separate reaction-plane detectors 
positioned in the range 1.0 < \rj\ < 3.9 show good agreement, indicating the absence of signiffcant 
Ari-dependent nonflow correlations. Sizable values for V4(pt) are observed with a ratio "jS i^^'"!^,'""*? ~ 
0.8 for 50 < A'part ^ 200, which is compatible with the combined effects of a finite viscosity and 
initial eccentricity fiuctuations. For A'^part ^ 200 this ratio increases up to 1.7 in the most central 
collisions. 

PACS numbers: 25.75.Dw, 25.75.Ld 



The discovery of large azimuthal anisotropy at the Rel- 
ativistic Heavy Ion Collider (RHIC) is a key piece of 
evidence for the creation of dense partonic matter in 
ultra relativistic nucleus-nucleus collisions [H, [1]. With 
sufficiently strong interactions, the medium in the colli- 
sion zone can be expected to locally equilibrate and ex- 
hibit hydrodynamically driven flow 3-5]. The momen- 
tum anisotropy results from an initial "almond-shaped" 
collision zone produced in noncentral collisions 0, |j| • It 
is now routinely characterized, at midrapidity, by the 
even order Fourier coefficients Vn = <^e*"('^p~*"P^) ,n = 
2,4,..., where 0p is the azimuthal angle of an emitted 
particle, $rp is the azimuth of the reaction plane and 
the brackets denote averaging over particles and events. 

At the highest RHIC collision energy of 
y^SNN = 200 GeV, differential elliptic flow measurements 
V2{pt) (for transverse momentum < 2.5 GeV/c) and 
?;2(iVpart) have been measured for a broad range of 
centralities or number of participants A^part- These 
data are found to be in accord with calculations that 
model an essentially locally equilibrated quark gluon 
plasma (QGP) having little or no viscosity d, 
Quark number scaling of elliptic flow data (suggestive 
of partonic degrees of freedom in the collision zone) is 
observed for a broad range of particle species, collision 
centralities and transverse kinetic energy [l3| • Small 
violations of the scaling of t'2(.^part) with the initial 
eccentricity of the collision zone e, suggest a strongly- 
coupled low- viscosity plasma (47r^ ~ 1 — 2 for the ratio 
of viscosity rj to entropy density s) in energetic Au+Au 
collisions 11-[T3|. Nonetheless, the degree to which the 
QGP is thermaliz ed |14i . and whether it is strongly or 
weakly coupled [1, Il5l|. is still being debated. 

Recent theoretical studies indicate that the hexade- 
capole flow harmonic V4 is a more sensitive constraint 
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on the magnitude of ^ and the freeze-out dynamics JL6.] , 
and the ratio -jj^ can indicate whether full local equi- 
librium is achieved in the QGP ,,17j. The role of fluc- 
tuations and so-called "nonflow" correlations is impor- 
tant for such measurements. It is well established that 
initial eccentricity fluctuations significantly influence the 
magnitudes of ^2,4 [18 22]. However, the precise role of 
nonflow, which leads to a systematic error in the determi- 
nation of 1)2.4 is less clear. Non-flow correlations among 
produced particles may arise from jets, whose influence 
is found to vary with pseudorapidity r] and pt [23l ] . This 
provides a tool to evaluate how jets influence the mea- 
surements presented here. 

We report precise measurements of charged hadron V2 
and V4 in Au-|-Au collisions at y^s^vjv = 200 GeV. The 
measurements were performed in the two PHENIX cen- 
tral arms (jr/j < 0.35) with respect to event planes ob- 
tained from four separate reaction-plane detectors in the 
range 1.0 < \ri\ < 3.9. Multiple event planes allow a 
search for possible A?7-dependent nonflow contributions 
that would influence the magnitude of 'U2,4, which may 
be crucial for reliable extraction of transport coefficients. 

The results reported here are derived from ~ 3.6 x 
10^ minimum-bias Au-|-Au events obtained at y^SNN = 
200 GeV with the PHENIX detector jH] during the 2007 
running period. The event centrality was determined via 
cuts on the analog response of the Beam-Beam Coun- 
ters (BBC). For each centrality selection, the number of 
participant nucleons A^part, was estimated via a Glauber 
model Monte-Carlo simulation [2^. The drift cham- 
bers and two layers of multi-wire proportional cham- 
bers with pad readout (PCI and PC3) were used for 
charged particle tracking and momentum reconstruction 
with azimuthal coverage Aip — tt/2 in the central region 
{\ri\ < 0.35). Tracks were required to have E/px > 0.1 
and a confirmation hit within a 2a matching window in 
PC3 and the Electromagnetic Calorimeters (EMCal) (E 
denotes the energy deposited in the EMCal). This mini- 
mized albedo, conversions and weak decay products. 

The event-plane method was used to correlate 
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FIG. 1: (color online) Event-plane resolution factors vs. TVpart 
for V2 (a) and 174 (b) measurements for the indicated event 
planes. 



the azimuthal angles 0p of the charged tracks in the 
PHENIX central arms (|?7| < 0.35) with the azimuth of 
the estimated second order event plane $2, determined 
via hits in the two BBCs and Muon Piston Calorime- 
ters (MPCs), and the two inner (i), outer (o) and com- 
bined (io) rings of newly installed Reaction-Plane Detec- 
tors (RXN). The two RXNs are situated at |z| =38- 
40 cm of the nominal crossing point and their inner and 
outer rings are comprised of twelve plastic scintillators 
(A0 = 7r/6 for each). The MPCs are PbW04 based elec- 
tromagnetic calorimeters with 27r azimuthal acceptance. 
The respective 77 coverage for these event-plane detec- 
tor pairs are 3.1 < |?7bbcI < 3.9, 3.1 < |77„pcl ^ 3.7, 
1.5 < |??RXNi| < 2.8 and 1.0 < |?7rxNo| < 1.5. For a given 
pair the detector, which is located at positive (negative) 
ry, is designated North (N) (South (S)). 

Charge-averaged values for the second and fourth flow 
harmonics were evaluated separately for each estimated 
event plane i as: 



(cos(2fc(0p - $1))) 



Res(^' 



fc = 1,2 



(1) 



2k! 



where the denominator represents a resolution factor that 
corrects for the difference between the true azimuth <i>Rp 
and the 2"'^ order estimate of the event plane. This 
estimate was obtained from the combined sub-events 
(North and South) for each detector pair. Resolution 
factors were evaluated via the three-sub-events method 



cos(2fc($^ - $^)))(cos(2fc($| - $™))) 



(cos(2fc(<i>^ - <i>™))) 



.(2) 



where j, I and m indicate event and subevent planes with 
disparate r/ values (eg., i = RXNio, / = MPCjv, and m — 
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FIG. 2: (color online) Comparison of 112 vs. A'part (a) and va 
vs. A'^part (b) for charged hadrons obtained with several re- 
action plane detectors for the pr selections indicated. Ratios 
for the pt range 1.2-1.6 GeV/c are shown in (c) and (d); the 
curves indicate ±10% and ±20% systematic error bands. 



BBCs). An advantage of this procedure is that, for any 
given centrality, it allows several independent estimates 
of Res(5'2.4) for each event plane. In turn, such estimates 
allow an evaluation of the systematic errors for Res(4'2 4). 
It is noteworthy that estimates for these correction fac- 
tors were also obtained (for k = 1 and 2) via the two- 
sub-events method [2^ [23|, which is regularly used for 
elliptic flow analysis. For RXN the difference between 
both methods is small for V2 i.e, ^ 1% for mid-central 
collisions and ^ 5% for the most central and peripheral 
collisions. For U4, it is ~ 2% for mid-central collisions 
and grows to ^ 7% and 20% in the most peripheral and 
central collisions respectively. 

Figure [1] shows the centrality dependence of (Res(\E'2)) 
and (Res(^'4)) for several event planes. Similar maxima 



are observed for TV, 



part 



200 with a falloff at lower and 
higher A^part- Measurements with the RXNio event plane 
benefit from about a factor of two (five) improvement 
in the resolution for V2 {vi) compared to prio r PHENIX 
measurements with the BBC event plane [26|. 

The systematic errors associated with the RXNjo res- 
olution factors for V2 {vi) are estimated to be less than 
2% (6%) for mid-central collisions but increase to about 
3% (10%) in the most central and peripheral collisions. 
Similar estimates were obtained for the RXNj and RXNo 
event planes. On average, those for the BBC and the 
MPC event planes are about a factor of two larger. Other 
sources, such as track cuts, are estimated to range from 
- 1 - 2% (3-4%) for PT > 0.5 GeV/c to -5% (10%) for 
the lowest pT values. 

Figures (Ha) and (b) compare the double differential 
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FIG. 3: (color online) pr dependence of V2 (a) and V4 (b) for 
charged hadrons for several centrality selections as indicated. 
The error bars only indicate statistical errors. 



flow coefScients W2,4(pt, -^part) for event-plane detectors 
spanning the range 1.0 < |?7| < 3.9. Within systematic 
errors, they agree to better than ~ 5% (10%) for V2 {V4) 
in mid-central collisions and approximately 10% (20%) 
in central and peripheral events (cf., ratios in Figs. [2][c) 
and (d)) independent of pr- This agreement indicates 
a reliable measurement free of significant Arj- and px- 
dependent nonflow contributions (for 3 GeV/c), 

which would affect V2 and W4 (very little influence is ex- 
pected from a possible Ary-independent long-range cor- 
relation [1^). Non-flow correlations, such as from dijets, 
would lead to a difference in the V2 (v^) values obtained 
with event planes determined at different rapidity gaps 
(A77) with respect to the central arms [2^. In the fol- 
lowing we utilize the RXNio event plane due to its good 
resolution. The associated systematic error for V2 (vi) is 
estimated to be w 3% (8%) for mid-central collisions and 
increase to about 7% (15%) in the most peripheral and 
central collisions. 

Figures [3] and HI summarize the results for elliptic and 
hexadecapole flow. The magnitude and trends agree well 
with those reported earlier [l|, Q . However they now ben- 
efit from a factor of five increase in statistics, as well as 
improved precision (~ 2) in the event plane. FigurestSJa) 
and (b) compare the measured charged hadron differen- 
tial V2{pt) and v^Ipt), as a function of centrality. In con- 
trast to the approximately linear dependence observed in 
Fig. ^a.) for pT < 1.5 GeV/c, the V4 data exhibit a non 
linear dependence on px compatible with the prediction 
from hydrodynamics that V4 oc z;| (29| . The large increase 
(~ x6) from central to peripheral collisions, refiects the 
expected increase due to the change in initial eccentricity 
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FIG. 4: (color online) 112 vs. A^part (a) and V4 vs. A''part 
(b) for charged hadrons for several pr selections as indicated. 
Panel (c) shows the ratio V4/{v2)'^ vs. A'part for the same pr 
selections. The open boxes indicate systematic errors for the 
selection 1.6 < pr < 2.4 GeV/c. The curves show calculated 
results for ideal hydrodynamics (I-Hd), I-Hd -I- eccentricity 
fluctuations (e-F), I-Hd -I- e-F -I- viscosity (77), and I-Hd + rj 
+ Gaussian Fluctations (GF) (see [2^). 



from central to peripheral events [13, llOl • 

Figured compares the ?;2(-^Vpart) (a) and W4(A'^part) (b) 
for several px selections as indicated. The A^part values 
are mean values evaluated for the centrality selections in- 
dicated in Fig. [31 Here, the data trends in (a) and (b) are 
strikingly similar albeit with a much smaller magnitude 
in (b) . The magnitude and trends with px and iVpart in 
FigS-HI^a) and (b) follow expectations for a Iwdrodynam- 
ically expanding low viscosity fluid 0, 0, @, illl - fisj . 

The ratio is shown as a function of A^part in 

Fig-HUc) for the same px selections used in (a) and (b); 
systematic errors are 4 — 5% for mid-central collisions 
and increase to 8-10% for central and peripheral colli- 
sions. Within errors, these data indicate that the mag- 
nitude of is essentially independent of px for the 
range 0.5 — 3.6 GeV/c, i.e. extending beyond the max- 
ima in Fig. [2I^a). An approximately constant ratio of 
value "^(P^'^p-*) ~ 0.8 is observed for 50 < iVpart < 200, 

t)^ (pT, Afpart) ^ pail , 

which is larger than the ratio ~ 0.5 for ideal hydrody- 
namics in the model of [23 | . The inclusion of eccentricity 



6 



fluctuations in this model, cause this ratio to exceed 0.5 
as shown by the dashed curve (from in Fig. |H^c). 
Viscosity from the hadron gas phase, in addition to a 
small value in the quark gluon plasma {Att- ~ 2) p^ . 
results in a further increase of this ratio as indicated by 
the dashed-dot curve p2l |. 

n„r "4(pT,Afp 
^Ul —7? 77 V 



ratio is smaller than the centrality- 



averaged value of 1.2 reported by STAR [3lJ. Part of this 
difference can be understood by averaging over our mea- 
sured centrality range (0-60%) yielding the value ~ 1.0. 
Comparison to STAR results [22] shows a 10% discrep- 
ancy for mid-central collisions, possibly reflecting differ- 
ences in the methods used to estimate Res(^'4). 

In more central collisions where iVpart ^ 200, ^ in- 

"2 

creases rapidly. Adding eccentricity fluctuations to ideal 
hydrodynamics causes a similar trend, indicated by the 
dashed curve in Fig.Hfc). Central collisions are the most 
sensitive because the eccentricity decreases as the overlap 
region becomes more symmetric. In order to reproduce 
the central data, the authors of [l^l introduced additional 
fluctuations shown as the solid line in Fig. |H[c), though 
the source of these fluctuations is as yet unspecified. 

In summary, we have presented differential measure- 
ments of W4 and V2 for charged hadrons obtained with 
four reaction-plane detectors at different A77 with re- 



spect to the PHENIX central arms. There are no sig- 
nificant Arj- and pr-dependent nonfiow contributions for 
Pt ^ 3 GeV/c in the centrality ranges of our study. Con- 
sequently there are no significant systematic errors from 
jets on the event-plane determinations or values of V2 
and V4. The ratio ~ 0.8 for 50 < TVpart < 200 

is essentially independent of pT, consistent with the ef- 
fects of finite viscosity and eccentricity fiuctuations. For 
A^part ^ 200 the ratio increases up to 1.7 in the most 
central collisions. The precision of these data provide 
stringent constraints for further theoretical modeling and 
more detailed extractions of the transport properties of 
hot and dense partonic matter. 
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